Enlarging and cooling the Neel state in an optical lattice 
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We propose an experimental scheme to favor both the realization and the detection of the Neel 
state in a two-component gas of ultracold fermions in a three-dimensional simple-cubic optical lattice. 
By adding three compensating Gaussian laser beams to the standard three pairs of retroreflected 
lattice beams, and adjusting the relative waists and intensities of the beams, one can significantly 
enhance the size of the Neel state in the trap, thus increasing the signal of optical Bragg scattering. 
Furthermore, the additional beams provide for adjustment of the local chemical potential and the 
possibility to evaporatively cool the gas while in the lattice. Our proposals are also relevant to 
efforts to realize other many-body quantum phases in optical lattices. 

PACS numbers: PACS numbers: 71.10.Ca, 71.10.Fd, 71.10.Hf,37.10.Jk 



A. Introduction 

Cold atom experiments provide a uniquely versatile 
platform for realizing and probing strongly correlated 
quantum phases of matter. However, no experiment to 
date has measured a phase in an optical lattice whose 
ordering is set by a magnetic scale such as the superex- 
change energy. Experiments have realized Mott insu- 
lators of both bosons [T] and fermions [H [3], but the 
temperatures achieved are higher than those required for 
magnetic ordering |3HZ1- Furthermore, the lack of a heat 
bath imposes restrictions on experimental schemes |8l- 
lllj . in which the entropy must be pushed out from the 
center of the trap where the phase of interest is realized. 
In the case of gapped phases, equilibration is impeded by 
the long timescales for heat and mass transport |T2l [E] . 
Finally, an experimental setup should strive to have the 
phase of interest occupy as large a region of the trap as 
possible to enhance the detectability of the ordering, for 
example, by Bragg scattering of light [Uj . 

In this work, we propose an all-optical scheme that 
addresses these issues for the Neel phase of ultracold 
fermions in a simple-cubic optical lattice, and discuss its 
relevance to other efforts to realize strongly-correlated 
many-body quantum phases. The objective is to maxi- 
mize the size of the phase of interest in the trap to en- 
hance the Bragg signal, realize the Neel phase in the re- 
gion of parameter space that was previously calculated to 
have maximal superexchange interactions |15) . and pro- 
vide a setup that will allow for cooling when the center 
of the trap becomes a Mott insulator, for which heat and 
mass transport are inhibited. We show that these ob- 
jectives can all be met simply by introducing three com- 
pensating laser beams on top of the three retroreflected 
lattice beams. These compensating beams have different 



Gaussian beam waists than the lattice beams and are 
oppositely detuned, so that if the lattice beams generate 
an attractive potential for the atoms, the compensating 
beams are repulsive, and vice-versa. The compensating 
beams allow the overall chemical potential of the system 
to stay in the gap of the phase of interest over a larger 
fraction of the cloud. We propose to do this in a dif- 
ferent manner than was analyzed in Rcf. [8^, where the 
trap potential was flattened by making it either quartic 
or something close to a square well. In both of these cases 
the walls of the trap are made steeper when the bottom 
is made flatter, and as a result the number of atoms in 
the outer gapless part of the cloud is reduced, making 
the system very sensitive to variations in the total atom 
number. In our setup, the confining potential is smooth 
and decays as a Gaussian at large distances. The trap is 
filled so that there remain many atoms in the outer gap- 
less "reservoir" parts of the cloud, as discussed below. 
Thus the system will not be sensitive to small variations 
in the total atom number. Furthermore, our setup allows 
for direct evaporative cooling of the system while in the 
lattice. 



B. The model 

We consider two-component ultracold fermions of mass 
TO in a simple-cubic optical lattice interacting repulsively 
via a Feshbach resonance, but far enough from the Fes- 
hbach resonance to apply the first order Born approx- 
imation. We call d the lattice spacing, the s-wave 
scattering length, and measure energy in units of the 
recoil energy Er = h'^kj^/{2m) = h'^n'^ /{2m,cP) where 
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FIG. 1: (Color online) Phases and potentials in the trap, (a 
Schematic representation of the distribution of phases. The 
three lattice beams run along the x-, y- and z-axes. They 
meet and form a simple-cubic lattice at the center. The anti- 
ferromagnetic (AF) Mott-insulating phase with filling of one 
atom per lattice site is found here. Surrounding the AF phase 
in the region where all three lattice beams are non-negligible 
is a paramagnetic Fermi gas phase Ri ("Reservoir I"). The 
regions where only one of three lattice beams has significant 
intensity are denoted Rn. The region where all lattice beams 
are negligible make up Rm. Atoms can be contained in Rm 
by an additional confining Vext{j')- |(b)| Total potential at 
z = 0: V{x,y,z — 0). Compensating beams, of optical po- 
tential opposite to the lattice beams, run along the x-, y- and 
z- axes to maximize the size of the AF region. 



kfj = TT /d is the recoil momentum. The Hamiltonian is 
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FIG. 2: (Color online) Plot of the Hartree estimate of the 
effective antifferomagnetic exchange Jh in a simple-cubic op- 
tical lattice with lattice depth Vq*"*, for interaction strengths 
aa = O.ld (blue) and as = 0.06d (red) where d is the lattice 
spacing. Jh is twice the difference in energy per bond between 
the antiferromagnetic and ferromagnetic states, obtained in 
the Hartree approximation (see [T5| for details). A large Jh 
leads to a large ordering temperature, and fast time scales for 
equilibration and entropy transport. The highlighted points 
(orange) correspond to the maximum of Jh as a function of Vb 
for a given interaction strength a^/d. Employing the local- 
density approximation, the trap shown in Fig. [l] is simple- 
cubic along the diagonal directions with lattice depth Vq*"® 
decreasing as one leaves the trap center. Thus by arranging 
to have VJJ^"® be close to the value that maximizes Jh over a 
large portion of the trap, the conditions for realizing the Neel 
phase are optimized. 



/3o-(r) = ^J.(r)5'(j(r) are the density operators; and V{r) 
is the total potential felt by the atoms. The first order 
Born approximation gives g = ATrh^as/m. 

The total potential V{y) is composed of an external 
potential Vexti^), three lattice beams, and three compen- 
sating laser beams. The external potential Vea;i(r) (which 
may be zero) is generally provided by optical dipole forces 
and varies over a length scale much larger than all the 
other length scales in the system. The three lattice beams 
and compensating beams are oriented along the x, y, and 
z axes. We call VLc{x;y,z) = VL{x;y,z) + Vc{x;y,z) 
the sum of an optical lattice along the x direction, pro- 
duced by a retroreflected Gaussian laser beam, and a 
non-retroreflected Gaussian beam along the x direction, 
that serves to partially compensate the overall average 
potential of the lattice beams: 

VL{x-y,z) = tV^l exp (- X ) ^i^' i^R^) (2) 



2(y2 _|_ ^2-) 

Vc{x\y:z) = ±Vbcexp( — ^ ) 



(3) 



where r = {x,?/, z}; ^cr(r) {^\{r)) is the fermionic an- 
nihilation (creation) operator of spin a at position r; 



The intensities Vol and Vqc of the lattice and compen- 
sating beams are positive. The upper signs correspond 
to having attractive lattice beams and repulsive compen- 
sating beams, the lower signs to the opposite situation. 
Generating a three-dimensional simple-cubic optical lat- 
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FIG. 3: (Color online) Phases along the diagonal direction (schematic). The horizontal axes denotes the radial coordinate 
r — \J -\- -\- . The red lines are schematic depictions of the potential along the diagonals: F(±a;, ±j/, ±2). The dashed 
lines indicate the value of the overall chemical potential \i. The red stars denote the lowest state of the system, within a local 
density approximation. If ^ lies below the red stars, the system is in the vacuum. The shaded blue regions correspond to 
the antiferromagnetic Mott insulator phase (AF). For above the red stars but outside the shaded blue region, the gas is a 
paramagnet, which we call Ki where the lattice potential is non-negligible. Rm is the region where all lattice potentials are 
negligible and is occupied only if /x > 0. There are four general classifications: (a) attractive lattice without compensation; 
(b) repulsive lattice without compensation; (c) attractive or repulsive lattice with compensation and Rm is non-empty; (d) 
attractive or repulsive lattice with compensation and Rm is empty. These classifications are depicted more quantitatively in 
Figures |4] to [S] With compensation, the AF phase has its largest possible size when /i is set so it coincides with the point 
where the Hubbard gap closes; this case is shown here. In the case of (d), ^ is chosen to be just below the zero of energy. This 
situation enables evaporative cooling, as particles with energies above zero correspond to excitations and are able to leave the 
trap. In (c), the system is in thermal contact with the outside reservoir Rm. The temperature of the system in equilibrium is 
set by the temperature of Rm, which itself may be evaporatively cooled during the experiment. 



tice requires three copies of the retroreflected lattice and 
compensating beams, in the three orthogonal directions. 
The lattice beams alone lead to a simple-cubic optical 
lattice in a region of space that is limited by the Gaus- 
sian profile of the beams, wl and wc are the waists of 
the lattice and compensating beams, respectively, and we 
define their ratio as a = wl/wc- The total potential is 
given by 



V{v) = V,Ar) + VLc[x\y,z) 



(4) 



The position and size of the different phases in this 
trap are determined by combining a calculation of the 
phase diagram of a homogeneous system with the lo- 
cal density approximation (LDA). In previous work |15| , 
two of us studied the phase diagram of this model in 
the Hartree approximation at zero temperature with a 

potential V^LiySin^^kjix) + sin^ikny) + siv^{kfiz) \ . To 



relate it to the present work, we assume that kj^wc S> 1 
and kf(WL 3> 1, so that the Gaussian envelopes of the 
potentials vary much more slowly than the lattice spac- 
ing. Under the LDA, the potential at each point in the 
trap is a sum of sinusoidal potentials plus an overall shift 
^■con{x,y, z) that is the local minimum of the lattice po- 
tential: 

V{v) = ^iconix, y, z) + Vo.j:{y, z)sin'^{k]ix + (j)^) 
+Voy{x, z)sin^{kiiy + (fly) + Vbz(x, ?/)sm^(fc_RZ -f 0^) (5) 

where (f)x,(t>y,4'z are relative phases which are unimpor- 
tant in the LDA. We neglect the spatial variations of the 
lattice amplitudes and ^con over the period of the lattice 
oscillation, but we do account for their variation over the 
longer scales given by the beam waists wl and wc- We 
choose the magnitudes of the sinusoidal parts of the po- 
tential to be positive: Vox{y,z), Voy{z,x), Voz{x,y) > 0. 
This can be done whether the lattice beams are repulsive 
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Blue-detuned lattice without compensation 

FIG. 4: (Color online) Phases without compensation and with 
aa — O.ld. (a) attractive, and (b) repulsive lattice beams. The 
horizontal axis is the local lattice depth VqI'^^ which depends 
monotonically on the distance r from the trap center along 
a {±1,±1,±1} direction. The zero of energy for the vertical 
axis corresponds to the potential in region Rm, where all the 
lattice and compensation lasers are negligible. The red stars 
denote the bottom of the band along a diagonal, and corre- 
spond to the lowest-energy single-atom Bloch state when the 
system is empty. The overall chemical potential jj, is denoted 
by the horizontal black dashed lines. If ^ < 0, Rm is empty; 

below the red stars corresponds to vacuum. The shaded 
blue regions correspond to the AF phase, /i outside of the 
blue region and above the red stars corresponds to a param- 
agnetic Fermi gas in reservoir Ri . The green plus signs denote 
the bottom of the band in reservoir Ru, where only one lat- 
tice beam is present. Ru is empty when n is below these. 
Since the laser potentials in the different directions are addi- 
tive, the band bottom along the x, y or z axis in Ru (green 
plus signs) is one third of the band bottom along the diagonals 
(red stars). The chemical potential is chosen so that the AF 
phase appears around the optimal lattice depth (Vb = AEr 
for fls = O.ld) calculated in the Hartree approximation [15) . 
The AF phase occupies a relatively narrow region of values 
of lattice depths, which translates into a narrow region of the 
trap. As the lattice gets deeper, the bottom of the band gets 
pulled down (up) for attractive (repulsive) lattice beams. 



or attractive by an appropriate choice of /iconix , y , z) , 
4>x, 4>y and 4>z, since for example — Vbs(j/, 2;)sm^(fc/jx) = 
-yoxly-, z) + {V()x{y, z))sin^{kiiX - 7r/2), and we can ab- 
sorb -VQx{y,z) into ficonix,y,z) and -7r/2 into 0^. 

The density in the trap is set by choosing an overall 
chemical potential /i for the system, assuming the system 
is at global equilibrium at zero temperature. Within the 
LDA, the system sees a potential Vox{y,z)sin'^{k]ix) + 
Voy{x, z)si'n?{kiiy) + Vqz[x, y)sin?{kjiz), and has a local 
chemical potential given hy fi — ^j,conix,y, z). At each 
point in the trap the Hartree calculation takes the local 
chemical potential and lattice intensity and returns lo- 
cal properties such as density, staggered magnetization 
and the Mott-Hubbard gap. The Hartree calculations 
were restricted to the case of equal lattice intensities, 
Vox = Voy — Voz, which occurs along 4 straight lines in 
the {1,±1,±1} spatial directions. We call Vjj*''^(r) the 
intensities of the lattice beams, and n^l^^ the chemical 
potential shift due to the lasers along the {±1,±1,±1} 
directions, where r — x'^ + y"^ + z-^ . They have the fol- 
lowing form: 

V,'rir)^VoLexpi~^J (6) 
Attractive lattice beams : (7) 
t^tn'ir) = i{V^cexp{-^) - VoLexp{-^)){%) 
Repulsive lattice beams : (9) 
Mcr(0 = 3( - V,cexp{-^)) (10) 

We choose the zero of energy so Vext = in the lat- 
tice (indeed we assume Vext is negligible in regions 
where other laser potentials are sizeable). We define 
a dimensionless parameter j3 to characterize the ratio 
of intensities of the lattice and compensating beams: 

{Vac/ER)exp[~^) = fi({yaL/ER)exp{-^^)) , so 
that 

Voc/En,^ P[VaL/ERr\ (11) 

Thus /3 gives the ratio of the intensities of the compensat- 
ing and lattice beams at the point where V^^^{r) — Er. 
The potential on the z = surface is depicted schemati- 
cally in Fig. [T(b)| 

C. The Neel state and its reservoirs 

The Hartree calculation in the lattice gives the regions 
of parameter space where the ground state is the Neel 
antiferromagnetic phase (AF), and the highest-energy 
occupied and lowest-energy unoccupied Hartree single- 
atom states provides an estimate of the Mott-Hubbard 
gap (neglecting spin- wave corrections). If the chemical 
potential lies anywhere within this gap, the phase will be 
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AF. Surrounding this AF phase is a "reservoir" of atoms 
in a paramagnetic Fermi gas. We want the atoms in this 
reservoir to be mobile so that they can carry away en- 
tropy from the part of the trap containing the AF phase. 
Thus, we want the optical lattice to remain relatively 
weak in the reservoir. 

We can distinguish between three different types of 
reservoir (Fig. [T]): Rj is in the three-dimensional part 
of the lattice, i.e. in the region where all lattice and 
compensating beams are non-negligible; Rn are the 6 
reservoirs in the regions where the beam intensities are 
appreciable in only one direction, corresponding to taking 
one of coordinates |a;|, \y\ or \z\ large compared to the 
beam waists while leaving the other two small enough to 
remain within the beam; and Rm is the reservoir outside 
of all of the beams. 

The parameters in the potential allow for significant 
freedom in tailoring the distribution of phases in the trap. 
As one goes along one of the diagonal directions away 
from the origin, the amplitude of the lattice decays ac- 
cording to Eq. [6j which sets the local effective amplitude 
^(faag ^j^g simple-cubic optical lattice. There is no need 
to specify the waist of the lattice beam within the LDA, 
as this simply sets the linear size of the different phases 
in the trap. The parameters that must be chosen are the 
ratio of beam waists a and the ratio of the intensities of 
the lattice and compensating beams /3, as defined in Eq. 
??. 

The chemical potential /i can be directly related to 
the density in the region Rm where all laser potentials, 
except possibly for Vexti'r), are zero (see Fig. T]). If ^ < 
V'ea;t(r) then Rjii is empty. For /i > Vexti^) the Hartree 
approximation gives the local density in Rm: 



TT d 



(12) 



D. Optimizing the parameters 



The parameters of the system should be chosen to max- 
imize the size of the phase of interest, and create optimal 
conditions for the realization of the phase. The lattice 
depth which maximizes the Hartree estimate of the ef- 
fective AF exchange Jh, thus giving the fastest equili- 
bration timescales and maximum Neel temperature, for 
a given interaction strength a^, was estimated in previ- 
ous work [TH], and is plotted in Fig. [2] as a function of 
lattice depth for the interaction strengths considered in 
this work: = O.ld and = 0.06d. We also found 
that the lattice depth which maximizes the entropy of 
the Neel state, in a calculation that neglects terms in the 
Hamiltonian beyond the Hubbard model [TBI, is close to 
this optimal AF exchange lattice depth. Therefore, the 
center of the trap should be at a lattice depth close to 
this optimal lattice depth. The AF phase should also 
occupy as large a volume as possible in the trap. 

Figure |3] schematically depicts the phases encountered 
in going along a diagonal from r = 0, where the lattice is 



deepest, to the edge of the trap, where the lattice depth 
goes to zero. Figures [3](a) and (b) show the cases of 
lattice potentials with no compensation, while Figs. |3][c) 
and (d) show two cases of compensated lattice potentials. 
Without compensation the chemical potential within the 
AF phase depends strongly on lattice depth, and conse- 
quently, the AF phase occupies a narrower region of the 
trap. One desired effect of compensation is to fiatten the 
chemical potential of the AF phase with varying lattice 
depth, in order to enlarge the AF region. In addition, as 
shown in Figs.jsjc) and (d), compensation allows adjust- 
ment of /i with respect to E = (the potential outside 
of the lattice). For fi slightly less than 0, Rm is empty, 
and atoms may evaporate from the edges of the lattice, 
while for fi slightly above 0, Rm is occupied and forms a 
thermal reservoir that may itself be evaporatively cooled 
and equilibrate with atoms at the lattice edge. 

More quantitatively, the results of our calculation of 
the phases are shown in Figures |4] to [8] In these plots, 
the horizontal axis is the lattice depth along a diagonal 
from the edge of the trap to the center. The center of 
the trap can be chosen to be anywhere along this axis 
(or beyond). Figure |4] shows the cases with no compen- 
sating beams. The smallness of the AF region due to 
the strong variation of the AF phase with lattice depth 
is readily apparent. To make the AF phase occupy a 
larger fraction of the trap, the compensating beam must 
shift the chemical potential of the AF phase so that it 
remains constant with varying lattice depth as shown in 
the following figures. 

Three different scenarios can be engineered for the fill- 
ings of the reservoirs, which we consider for as = O.ld 
in the following plots: Rm non-empty (Fig. [51; Rm 
empty (Fig. [6|; and Ru and Rm empty (Fig. 7^. To 
get a sense of the dependence on interaction strength, we 
plot the phases for — Q.OGd when all three reservoirs 
are occupied in Fig. [8| 

A qualitative understanding of why different waists for 
the lattice beams and the compensating beams is advan- 
tageous is illustrated by Fig. [9] Unequal beam waists 
will cause the strength of the compensating beams to 
grow as a power law in the strength of the lattice beams, 
as expressed by Eq. ??. For deep lattices, the bottom 
of the Mott gap becomes narrowly separated from the 
bottom of the lowest Bloch band. The objective of the 
compensating beams is to keep the chemical potential 
inside the Mott gap as the lattice depth varies. In the 
absence of the compensating beams, the dependence of 
the bottom of the Mott gap on lattice depth is reason- 
ably well described by a power law, so the compensating 
beams can flatten the Mott gap, thus achieving our ob- 
jective. Furthermore, we expect from the behavior of 
the band bottoms as a function of lattice depth shown 
in Fig. |9] that we need a > 1 for attractive, and a < 1 
for repulsive lattice beams to obtain as flat a Mott gap 
as possible with the given setup. The special case where 
the lowest order confinement is quartic [S] is achieved by 
setting Vol/Vqc = ct^- While this choice of parameters 
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FIG. 5: Phases for aa = Q.ld with compensating beams, for 
(a) attractive, and (b) repulsive lattice beams, with > 
chosen to give density is nd?' = 0.1 in Rm . The parameters 
are (a) /3 = 0.379, a = 1.13; and (b) /3 = 0.705, a = 0.81. 
The lines and symbols are defined as in Fig. [4] 



flattens the bottom of the potential, it does not maximize 
the volume of the Neel phase, since that is achieved by 
flattening the Mott gap instead, as we propose. 

The Hartree approximation we have used in this work 
is a mean-field approximation and the real system will be 
quantitatively somewhat different from these Hartree es- 
timates, due to nontrivial fluctuations and correlations. 
However, we do not require precise numerical results to 
show the effectiveness of the presented scheme. The 
Mott-Hubbard gaps are large enough in the region of 
interest that quantitative changes in the precise values 
will not destroy the general qualitative features that this 
scheme relies upon. The solid (blue) lines in Figures |4] 
to [8] bound the AF Mott insulating phase. The Hartree 
approximation probably overestimates the range of sta- 
bility of this phase, since it does not include all quan- 
tum fluctuations. But if the AF phase is quantitatively 
a little smaller in these figures than what we show, the 
chemical potential and compensating potential can still 
be adjusted to enlarge the AF phase to fill the entire cen- 



FIG. 6: (Color online) Empty Rm reservoir with as = O.ld. 
(a) Attractive lattice with /3 = 0.315, a — 1.17; (b) Repulsive 
lattice with /3 = 0.80, a = 0.77. 



tral region of the lattice, and allow continued evaporation 
from the AF phase within the lattice. 



E. Maximizing the Bragg signal 

The AF phase can be directly detected using Bragg 
scattering of near-resonant light [M]. The Bragg sig- 
nal from scattering off the up spin density, for example, 
is proportional to the volume of the AF phase and to 
the square of the Fourier transform of the up spin den- 
sity at momentum (kn, k;^, kf;). We have shown that 
the volume of the AF phase may be maximized by vary- 
ing the relative intensities and waists of the lattice and 
compensating beams. While cooling and equilibration 
times are minimized at relatively low lattice depths, the 
Bragg signal is enhanced at deeper lattices, for which 
quantum fluctuations due to the site being doubly- or 



unoccupied are weaker. Figure 10 shows a plot of the 
Fourier intensity in the ground state as a function of 
lattice depth for = 0.1c? and — O.OGd. The lat- 
tice depths which maximize AF superexchange in the 
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FIG. 7: (Color online) Both Rn and i?/// are empty, with 
as — O.ld. (a) Attractive lattice with /3 = 0.30, a — 1.16. In 
this case, the AF phase does not go all the way to 2.4i5_R. (b) 
Repulsive lattice with /3 = 1.0, a = 0.74. For this setup in 
both (a) and (b), the evaporating atoms go out the "beams" 
through the empty Rii^s and not through Rm, which is at 
too high an energy relative to fi. 



Hartree approximation [T5] are indicated. We see that the 
Bragg signal is maximized by going to deeper lattices and 
stronger interactions. Therefore, one must compromise 
between the conditions that minimize the time scales for 
equilibration and cooling and those that maximize the 
Bragg signal. 

One way to strengthen the Bragg signal is to cool and 
equilibrate at the relatively low lattice depths that max- 
imize superexchange, but before performing Bragg scat- 
tering ramp up the lattice depth [T71 |T5] at a speed that 
is sufficiently adiabatic to reduce the quantum fluctua- 
tions. Since these fluctuations arise from virtual pairs 
of empty and doubly-occupied sites due to the superex- 
change process, one should be able to remove them pro- 
vided that the lattice ramp is adiabatic with respect to 
the Mott-Hubbard gap. The spin-wave zero-point fluctu- 
ations that are present in the corresponding Hciscnbcrg 
model are not strongly reduced in the limit of a deep 
lattice, while the virtual vacancies and doubly occupied 
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FIG. 8; (Color online) Similar situation to Fig.js] with a lower 
scattering length: as — 0.06ci. (a) Attractive lattice with 
a = 1.13, /3 = 0.389; (b) Repulsive lattice with a = 0.79, 
/3 = 0.72. 



sites are strongly suppressed. At nonzero temperature 
thermal fluctuations will also produce real empty and 
doubly-occupied sites. The lattice ramp will bias the hop- 
ping of these thermally-excited site defects and possibly 
produce heating unless the compensating beams are care- 
fully ramped together with the lattice to eliminate such 
forces. The ramp should be fast enough to shut down the 
hopping in order to freeze in these thermal excitations be- 
fore they can recombine and release their energy, which 
becomes increasingly high compared to the energy of spin 
fluctuations as the lattice strengthens. This suggests an 
optimal ramp rate is the fastest possible while remaining 
adiabatic with respect to the Mott-Hubbard gap in the 
bulk of the AF phase. The ramp will be nonadiabatic 
near the outer edges of the AF phase, where the lattice 
is weak and the Mott gap is initially very small, causing 
some of the Bragg signal to be lost, but since the Neel 
ordering was initially very weak there the gain by en- 
hancing the Bragg signal over the bulk of the AF phase 
should outweigh this loss near its edges. The precise 
balance between these various nonequilibrium dynamical 
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FIG. 9: (Color online) Bottom of the lowest Bloch band 
along the diagonal directions, i.e. for a potential V{x, y, z) = 
Vq^^ {sin^ {knx) + sin^{kuy) + siii? {kaz)) for (a) attractive 
lattice beams {Vq^^ < 0) and (b) repulsive lattice beams 
(Vq*"* > 0). The energy difference between the bottom of 
the lowest Bloch band and the bottom of the Mott gap de- 
cays quickly with increasing lattice depth. Therefore, to keep 
the chemical potential in the Mott gap one has to flatten 
the bottom of the band, for deep lattices, (a) Full red line: 
bottom of the band in the case of attractive lattice beams. 
The dashed (dot-dashed) line gives the asymptotic behavior 
at deep (weak) lattice. For a weak lattice, first-order pertur- 
bation theory gives the bottom of the gap to be —3VqI'^^/2. 
For a deep lattice, using the harmonic and lowest order an- 
harmonic terms for the wells in the lattice, one gets that the 

band bottom goes like -SF*"'' + ^{^/v^l"^ /Er - 1/4)Er. 
The band bottom goes down superlinearly in lattice depth, 
so the compensating beam intensity must therefore grow su- 
perlinearly in VqI""^ . (b) Full blue line: bottom of the band in 
the case of repulsive lattice beams. The dashed (dot-dashed) 
line gives the asymptotic behavior at deep (weak) lattice. For 
a weak lattice, the band bottom goes as ~ ■SVqI"'^ /2, while for 

a deep lattice, it becomes ~ 3(y^V^*''^/i5fl — l/4)i5_R. There- 
fore, the band bottom grows sublinearly in lattice depth, so 
the compensating beam intensity must grow sublinearly in 

considerations is a challenge, and deserves further study. 

Bragg scattering relies on the antiferromagnetic order- 
ing being along the spin direction set by the "up" and 
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FIG. 10: (Color online) Square of the Fourier transform at 
momentum (kn, kn, kn) = (vr/d, vr/d, 7r/d) of the up spin den- 
sity, where d is the lattice spacing, and v is the volume of the 
system, for two interaction strengths: as — O.ld (blue) and 
as — 0.06d (red). The Bragg signal obtained from scatter- 
ing light off the up spins will be proportional to this quan- 
tity. At infinite lattice depth, it becomes 1/4, as the Fourier 
transform of infinitely localized particles on the FCC lattice 
at this momentum is 1/2. The highlighted (orange) points 
are at the corresponding lattice depth where the AF superex- 
change is maximized. While a deeper lattice and stronger 
interactions lead to more localized particles, and therefore a 
stronger Bragg signal, deeper lattices also lead to smaller su- 
perexchange and therefore smaller ordering temperatures and 
longer timescales for heat transport. Both objectives can be 
met by cooling with a weaker lattice and then ramping up the 
lattice before performing Bragg scattering. 



"down" hyperfine states, which we call the z-direction. 
The local up and down spin populations may not be 
precisely equal, however, because of fluctuations in the 
initial conditions. This local spin polarization along the 
z direction produces canted antiferromagnetism [191 120) , 
in which the AF order is tilted only slightly away from 
the xy plane. In this case, a 7r/2 pulse before the Bragg 
measurement will tip the AF order up to a plane con- 
taining the z-direction, making it detectable with Bragg 
scattering [T4] . 

If the lattice depth is ramped up before Bragg scat- 
tering, the 7r/2 pulse should be applied before the lattice 
ramp for the following reasons: At deep lattices, the spin- 
spin interaction is greatly reduced and thus the compo- 
nents of the spins pointing in the xy plane on different 
lattice sites may dephase with respect to one another due 
to thermal and quantum fluctuations. However, spins 
pointing in the z-direction will not dephase since they 
are eigenstates of the single-atom Hamiltonian in a sin- 
gle deep well. Therefore, a 7r/2 pulse before ramping the 
lattice depths will partially prevent dephasing of the an- 
tiferromagnetic correlations that can occur at large Vql, 
by increasing the amount of the antiferromagnetic order 
that is along the z-direction. 
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F. Conclusions 

We have proposed a setup to facilitate both reahzing 
and detecting the Ncel state of two-component fcrmions 
in a simple-cubic optical lattice. We found that the in- 
troduction of compensating beams with a different beam 
waist allows for a significant growth of the Neel phase 
in the trap, and control over the different reservoirs that 
this state is in contact with. The ability to grow the size 
of the Neel phase in this simple setup relies on the ob- 
servation that the chemical potential of the Neel phase 
has a dependence on the lattice depth which is well ap- 
proximated by a power law. Since this is likely to be 
the case for other phases of cold atoms in optical lattices 
we expect that the proposed setup will confer similar ad- 
vantages to other attempts at realizing and probing such 
phases. One of the main challenges is realizing a setup 
where the system is able to shed its entropy, even as a 
gapped phase is forming and inhibiting transport. Typ- 



ically, present experiments rely on precooling the atoms 
and then adiabatically loading them in to the lattice and 
forming the phase of interest. Our proposed approach 
is to instead continue evaporative cooling as the lattice 
is turned on, by maintaining the chemical potential at 
a level that allows the phase to stay in contact with a 
reservoir of mobile atoms that is evaporatively cooled. 
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